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METHOD OF OPTICAL CHARACTER I SAT I ON OF MATERIALS WITHOUT 

USING A PHYSICAL MODEL 



DESCRIPTION 

5 TECHNICAL FIELD 

The invention relates to a method of optical 
characterisation of materials. 

This method allows to characterise thin or thick 
layers of these materials, which are created on 
10 substrates . The physical dimensions, that this method 
allows to determine, are: 

- the thickness of a layer of material; 

- the refraction index of this material; and 

- the absorption coefficient of this material. 

15 

Optical characterisation of materials is used for 
the chemical analysis of these materials (notably the 
study of absorption bands, of densif ication properties 
and of oxidation properties) ,. in the fields of 

20 microelectronics, sensors, biology, medicine, or to 
analyse the thickness of deposits of these materials. 

We refer to document [1] for examples of 
applications which, as for the other documents cited 
later on, is mentioned at the end of this description. 

25 The characterisation of the optical properties of 

a material is also useful when the material is 
structured later on (to create for example etchings or 
surface roughness) and the optical diffraction 
properties of the obtained structure must be calculated 

3 0 (see document [2] ) . 
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It can already be noted that the invention is 
particularly useful when the physical law applied by 
the complex refraction index of the material to be 
characterised is, a priori, unknown. 

5 

STATE OF THE PRIOR ART 

Remember that the optical measurements can have a 
variety of natures. 

It can be ref lectometric measurements. In this 
10 case, the intensity reflection coefficient of a 
structure is measured over a spectrum (i.e. an 
interval) of wavelengths [A m , A M ] . 

The incident angle of the illuminating light may 
be not zero. The reflection coefficient can be measured 
15 for several incidence angles 0. Let R (0, A, p) be the 
ref lectometric spectrum, where p is the polarisation of 
the incident beam and A the wavelength of the latter. 

Generally, the angle 0 is zero and the 
polarisation p unknown. In the event where 0 is not 
20 zero, this polarisation p must be known. Generally, the 
latter is of (S) or (P) type. 

It can also be ellipsometric measurements. The 
measured dimensions are then the real and imaginary 
parts of the ratio between the polarisation reflection 
25 coefficient (P) and the polarisation reflection 
coefficient (S) . 

We generally mark p=| pi exp { j A) this complex ratio 
(with j 2 = - 1 ) . The dimensions commonly handled are I p| , 
that we mark tan (W) , and cos (A) , or combinations of 
30 both. 
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For example, the variables derived from a phase 
modulation ellipsometer are the following: 

I s =sin(2¥) sin(A) and I c =cos(2¥). 
As for a standard ellipsometer, it supplies the 
5 following variables: 

a=(tan 2 ?-l)/(tan 2 ?+l) and p=cosA ( l-a 2 ) m . 
For generality purposes, we mark the handled 
dimensions S 1 and S 2 . 

The spectrums Si, ie [1,2], are measured over a 
10 range of wavelengths [A m/ A M ] . The incidence angle can 
be of any value. Several spectrums can be measured at 
different incidence angles in order to obtain a fuller 
spectrum. We mark the ellipsometric spectrum 
s (8, A) ={ Si (9, A) , S 2 (9, A) }. 
15 Additionally, goniometric measurements (reflection 

coefficient as a function of the incident angle) can be 
added to the measurement used for characterisation, so 
as to determine the thickness of the various layers, 
for one or several wavelengths. These measurements 
20 alone are insufficient as we want to determine the 
complex refraction index over a spectral range from A m 
to A M . 

So as to simplify the presentation, we mark W a 
set of ref lectometric and/or ellipsometric spectrums 
25 (and possibly goniometric spectrums for a few 
wavelengths) . 

Without omitting any generality, we describe, in 
this description, the operating mode of the methods of 
the prior art and of the invention only for the case of 
30 a single thin layer of material, created on a known 
substrate . 



B 14206.3 PV 



4 



The thickness of this layer is marked e and the 
complex refraction index of the material with the 
wavelength X is marked n* (X) . 

In this regard, remember that the real 
5 (respectively imaginary) part of this complex 
refraction index is marked n(X) (respectively k(X)) and 
called "refraction index" (respectively extinction 
coefficient" ) . 

Additionally, we mark an error function 
10 Er \p< 2 >) - for example the average quadratic gap - 

between two spectrums and W {2) . 

For example, we can take, when there are 
ellipsometric spectrums over several angles 0i, ie 
{l..n}, and a ref lectometric spectrum: 

15 Er <¥ U) ,W (2) ) = 

-^T[i£> 1 (,) (.pi).8 I «(. 1 .x)r + ^ , >(. l .x).s 2 m (. l .x)r " 
4* ( %)-r (2 w y CD 

with (X)={Sx (1) (9 if X) ,S 2 (1) (9i,X) ,R (1) (X) } 

20 ^ (2) (X)={S! (2) (6i,X) ,S 2 (2) (9 if X) ,R (2) (X) } 

and ie{l...n} 

Weighing factors can be introduced into the 
integral so that the error function can take into 
account the variations in the accuracy of - the 
25 measurements of the spectrums. 

The optical characterisation of layers of material 
is generally based around two applications. 
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The first application is the dimension control of 
the deposit of thin layers that are used in 
microelectronics . 

Usually we know the deposited material, i.e. we 
5 know the complex refraction index of this material at 
the first wavelengths used for the characterisation. 

The laws applied by the complex refraction index 
are either tabulated or approximated by known physical 
laws such as, for example, the Cauchy model, the 
10 Sellmeier model (see document [3] ) , the Forouhi laws 
(see document [4] ) , and the harmonic oscillator laws 
(see document [5] ) . These laws are defined by a finite 
number of parameters . 

For example, a Cauchy type law without absorption, 
15 with two parameters, is defined in the following 
manner : 

Re[n*(^)] = n(^)=a 0+ ^- 
Im[«*(Jl)]=k(A,)=0 

When we are certain of the value of the 
20 coefficients ai(i€{o,l}) but we do not know the 
thickness, a search algorithm is used in order to find 
the thickness which minimises the error between the 

measurement ¥ and the theoretical result ¥ taking into 

consideration the modelled index. 
25 The search algorithm can be, for example, the 

Simplex method, the Tabou search, the Levendt -Marquart 

method or the simulated annealing method (see chapter 

10 in document [6] ) . 

When the refraction index is approximate, the 
30 coefficients a± are integrated into the adjustment 
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process of W and W . The search for the coefficients ai 
constitutes a method of characterising the refraction 
index . 

However, when the law applied by this refraction 
5 index is unknown (it happens that the material is 
unknown or that it is poorly described by known 
physical laws) , this method remains approximate and the 
thickness risks being erroneous . 

The second application is the characterisation of 
10 materials. 

The employed method remains unchanged, except that 
the material is not really known. It is precisely the 
complex refraction index function that is the closest 
to reality that is targeted, 
is The type of law can be chosen via analogy with 

other materials. However, the law applied by the 
complex refraction index can be complicated, which is 
for example the case of a harmonic oscillator law. 



20 



[n(E) + jk(E)] 2 =l + l i 



i=lE + E- + jr. E-E- + jr. 

1 J 1 ^ I 



In the above expression, j 2 =-l and the refraction 
index and the extinction coefficient are expressed as a 
function of E and not of A, with E=1240/A(A in nm) . 
25 In this case, the coefficients of the oscillators 

are difficult to find if we do not know their value 
range. The search is difficult to automate, as the 
search algorithms can give erroneous results and the 
time wasted can be substantial. 
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There is an alternative to the search for 
coefficients: the point to point method (PAP). This PAP 
method proposes not to choose a physical law and to 
search the complex refraction index of the material for 
5 each wavelength X±, where i<= [l...n] , with X 1 =X m and 
A n =A M . 

For each Ai, a search algorithm tries to find the 
thickness, the index n(Ai) and the extinction 
coefficient k(Ai) which minimise the error between the 
10 measured values W (X±) and the theoretical result 

W (Ai,n(Ai) ,k(A ± ) ,e) . 

Such a method poses a problem because the 
different points ( X±, e , n ( X±) , k ( X±) ) are not necessarily- 
physical ly compatible among themselves: for example, 
is the found thickness can vary depending on the 
wavelength and the law applied by the complex 
refraction index, more simply called the index law, can 
have discontinuities . 

This method is usually only valid when the 
20 thickness is fully known and that the measurements are 
of good quality. 

PRESENTATION OF THE INVENTION 

25 The purpose of the invention is to resolve the 

previous inconveniences . 

The method which is the object of the invention 

allows to characterise a material without using a 

physical model, i.e. without using a physical law 
30 applied by the complex refraction index of the material 

under scrutiny. It is therefore especially useful when 

such a law is unknown. 
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This method constitutes an alternative to the 
aforementioned known methods of characterisation. It 
can be called "method of nodes" as it uses "nodes" , 
i.e. points having coordinates (Ai, n*i) , where n*i is 
5 the value taken by the complex refraction index with 
the wavelength X± and i. takes a limited number of 
values (whole numbers) . 

Precisely, the purpose of the invention is a 
method for optical characterisation of at least one 

10 layer of material in an interval A of values taken by a 
function a of an optical wavelength A, when A varies in 
an interval of wavelengths, this layer being created on 
a substrate, this method being characterised in that it 
comprises the following stages: 

15 1) we carry out a group of ref lectometry and/or 

ellipsometry measurements over the interval A, this set 
of measurements leading to a measured spectrum, marked 
and we choose the methods for calculating associated 
with the nature of the measurements and with the type 

20 of layer to be characterised; 

2) we choose m initial values ai ...a m of the 
function a, belonging to this interval A, m being a 
whole number at least equal to 1, and we define an 
interval B as being the set of points a of the interval 

25 ranging from the smallest to the biggest number among 
a x ...a m , when m is greater than 1, and as being the 
interval A when m equals 1; 

3) we choose m complex initial values of a complex 
refraction index n*=n+jk for the m points a±, i ranging 

30 from 1 to m; 
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4) when m is not 1, we choose an interpolation law 
which allows to calculate the refraction index n(a) of 
the material over the interval B, from the points (ot±, 
n±) , with n i =n(cxi) / i ranging from 1 to m, and when m 

5 equals 1, n(a) is taken equal to the number ni(ai) over 
the entire interval B; 

5) we choose M variable parameters, M being less 
than or equal to 2m+l; 

6) we choose an error function Er (W, ¥ ) which 
io characterises the difference between a measured 

spectrum W and a theoretical spectrum ¥ ; 

7) using a minimising function of Er {W, W ) with M 
parameters, we perform the following series of stages: 

a) by applying the interpolation law of (a^ni) 
15 over the interval B, we deduce n(a), a belonging to B; 

b) by using n(a) and the thickness e of the layer, 
and methods for calculating spectrums, we calculate a 

theoretical spectrum W (n(a),£); 

c) we compare W and T by using Er {W, W ) and, if 

20 Er (W, W ) is sufficiently small, i.e. less than a 
predetermined value e, or is minimal, we go to stage 
e) , otherwise we go to stage d) ; 

d) we make the M variable parameters vary so as to 

tend to the minimum of Er (W , ¥ ) , and we return to 
25 stage a) ; 

e) if Er W ) is less than e, we then obtain a 
set of M variable parameters, for which 

Er(¥, W (n(a,M),e)) is minimal and the refraction 
index is then taken equal to the last one obtained, and 
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if Er (W, W ) is greater or equal to e we go to stage 
8) . 

8) we increase the number m of initial values of 
the function a and we return to stage 2) . 
5 It is therefore, for example, possible to perform 

optical characterisation: 

— over an interval of wavelengths A, in this 
case the interval [Amin, Amax] , 

- or over an interval of inverse wavelengths 
io 1/A, in this case over an interval [(1/A)min, 

(1/A)max], where (1/A)min is equal to 1/ (Amax) and 
(1/A)max is equal to l/(Amin); 

- or over an interval of energies E (with E = 
hv = hc/A where h is the Planck's constant, c the speed 

15 of the light in the vacuum and v the frequency 
corresponding to (A) , in this case over an interval 
[Emin, Emax] , where Emin is equal to he/ (Amax) and Emax 
is equal to he/ (Amin) ; 

— or, more generally, over an interval [amin, 
20 amax] of values taken by a function a of the variable 

A. 

Furthermore, it is appropriate to note that the 
invention can be used to characterise a spectrum or a 
part of a spectrum. 
25 Each interpolation law can be chosen from among 

the linear interpolation laws, the cubic interpolation 
laws, the polynomial interpolation laws and the 
interpolation laws for example of "spline" function 
type. 

30 According to a preferred embodiment of the method 

which is the object of the invention, the spectrum is 
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evenly sampled in oc (A) , meaning that the initial 
values of the function a (see the aforementioned stage 
2) are evenly distributed over the interval A, the 
distribution of the nodes thus being homogenous. 
5 As previously seen, a (A) can be chosen from among 

A, 1/A and hc/A or any other function of A, where h is 
the Plank's constant and c is the speed of the light in 
the vacuum. 

Preferably, in the aforementioned stage 6) , we 
10 measure the error over an interest interval C which is 
included in the interval B or equal to this interval B. 

The M variable parameters can be the real parts of 
the refraction indexes at points ai, i ranging from 1 
to m, or the imaginary parts of these refraction 
15 indexes, or these M variable parameters can be 
constituted by the thickness of the material for which 
we are searching the refraction index. 

Another purpose of the invention is another method 
for optical characterisation of at least one layer of a 
20 material in an interval of wavelengths [A min, A max] , 
this layer being created on a substrate, this other 
method being characterised in that: 

- we carry out a set of ref lectometry and/or 
ellipsometry measurements, this set of measurements 

25 leading to a measured spectrum, marked 

— we choose m initial wavelengths Ai ...A m 
belonging to this interval, m being a whole number at 
least equal to 1, we associate a refraction index to 
each wavelength; 
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- we choose an interpolation law at least for 
the refraction index of the material, for wavelengths 
lying between the initial wavelengths A x ...A m ; 

- we choose M initial parameters, M being at 
least equal to m, namaly an initial refraction index ni 
for each initial wavelength Ai, 1 < i < m, the initial 
wavelengths being chosen so as to determine via 
interpolation at least the refraction index for any 
wavelength within the interval [A min, A max] , the 
couples (Ai, ni) being called nodes; 

- we choose ref lectometry and ellipsometry 
methods of calculation; 

- we also choose an error function Er, 
representative of the difference between two spectrums 
^ a and W 2 , the spectrums and W 2 being calculated or 
measured over a number of points greater than the 
number m of nodes; 

- using the m initial wavelengths, the M 
initial parameters and the interpolation law, we 

implement the following optimisation process: 

t. 

- we determine a theoretical spectrum, marked 

¥ , depending on the chosen methods of calculation, and 
on the index deduced via interpolation of its value at 
Ai, i ranging from 1 to m, over the spectrum [Amin, 
Amax] ; 

- we determine the error Er , W ) , between the 
measured spectrum and the theoretical spectrum; 

- we minimise this error by varying the 
position of the values of the unknown indexes and/or 
the thickness of the layer and/or the values of the 
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refraction indexes with initial wavelengths, and we 
obtain a spectrum; 

- we add other wavelengths to the initial 
wavelengths Ai ...Am, the added wavelengths constituting 

5 new nodes; 

- we repeat the method by choosing a number m' 
of initial wavelengths, m' being greater than m, and M' 
initial parameters, M' being greater than M, until the 
accuracy of each spectrum thus best represented is 

10 equal to a predetermined accuracy. 

In this case, according to a first specific 
embodiment, m is at least equal to 2; according to a 
second specific embodiment, m is equal to 1 and we 
choose equal initial refraction indexes. 

15 Also in this case, according to a specific 

embodiment, the material is non absorbent and the 
number M is equal to m, the extinction coefficient of 
the material being set equal to 0; according to another 
specific embodiment, M is at least equal to 2m, we 

20 furthermore choose an interpolation law for the 
extinction coefficient of the material, for each 
initial wavelength Ai, 1 < i < m, furthermore we choose 
an initial extinction coefficient ki, the initial 
wavelengths furthermore being chosen so as to be able 

25 to determine via interpolation the extinction 
coefficient for any wavelength of the interval [A min, 
A max] , and within the optimisation process, we 
minimise the error by also varying the values of the 
extinction coefficients at the initial wavelengths and 

3 0 the added wavelengths are furthermore placed so as to 
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best represent the spectrum of the extinction 
coefficient of the material. 

In the case of this other specific embodiment, m 
can be equal to 1 and we can choose equal initial 
5 refraction indexes and equal initial extinction 
coefficients . 

Still in the case of this other method which is an 
object of the invention, the layer of material can be 
thin, i.e. with a thickness less than the coherence 
10 length of the light used for measuring, we can choose 
an additional initial parameter, namely an initial 
layer thickness, and in the optimisation process we can 
minimise the error by also varying the value of the 
layer thickness; in an alternative, the layer of 
15 material can be thick, i.e. not thin, and M can be at 
most equal to 2 m; in another alternative, the 
thickness of the layer of material can be known with 
sufficient accuracy and M is at most equal to 2 m. 

The distribution of the nodes can be homogenous. 

20 

BRIEF DESCRIPTION OF THE DRAWINGS 

The invention will be better understood upon 
reading the description of the below embodiments, given 
25 by way of non-restrictive examples, making reference to 
the annexed drawings among which: 

- figure 1 is a diagrammatic view of devices 
allowing to characterise a layer according to the 
invention; 

30 — figure 2 shows the variations in the 

refraction index as a function of the wavelength, for a 
material obeying to a Cauchy law (curve I) and for a 



15 



material characterised according to the invention 
(curve II) ; 

- figure 3A (respectively 3B) shows the 
variations in the refraction index (respectively in the 
extinction coefficient) as a function of the 
wavelength, for a material obeying to a law with two 
harmonic oscillators (curve I) and for a material 
characterised according to the invention; and 

- figure 4 diagrammatically illustrates the 
parameters used in a generalisation of examples of the 
invention. 

DETAILED DESCRIPTION OF PARTICULAR EMBODIMENTS 

The invention proposes an alternative to the 
aforementioned standard methods. It allows to combine the 
coherence of a layer model (corresponding to a continuous 
index law and to constant physical thicknesses) , with the 
generality concerning the index law to be found (as in 
the PAP method) . Furthermore, the resolution is only 
limited by the resolution of the measured spectrum. 

In the method according to the invention, the 
spectrum of index n* (A) is characterised by: 

- a limited number of "nodes" which are points 
having co-ordinates (X if n±, k ± ) or (Ai, n*i) , with n ± =n 
(Ai) , ki=k (Ai) and n*i=ni+jki, where j 2 = -l, and 

- an interpolation law between the nodes, which 
can be, for example, linear, cubic, of "spline" type or 
polynomial (of any given degree) . 

This interpolation law allows to calculate, from the 
nodes, the refraction indexes and the extinction 
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coefficients for the wavelengths located between the 
nodes . 

For example, when the refraction index is 
characterised by a set of value for wavelengths Ai ...A m , we 
5 can use a linear interpolation between two wavelengths X± 
and Ai +1 to calculate the index n at the wavelength A (see 
document [6] chapter 3) : 

10 with Ai<A<A i+ i 

We can do the same for the extinction coefficient. 
When the number of nodes makes it possible, more 
complex interpolation formulae, bringing into play nearby 
nodes, can be used (see document [6] chapter 3) . 
15 A layer model is therefore characterised by a 

thickness e and a family of nodes. 

Hereafter we present an example of the method 
according to the invention. In this example, the 
measurements ¥ are constituted of a ref lectometry 
20 measurement R(A) and an ellipsometry measurement 
Si ,2 (9, A) where G is the incidence angle of the light 
beam that is transmitted onto the layer to be studied 
during measuring via ellipsometry. 

This layer is a thin layer so that the thickness of 
25 this layer is also a variable of the problem. Moreover, 
we suppose that only one layer is unknown, this layer 
being created on a known substrate. 

We will first of all briefly explain this example 
which uses an algorithm (algorithm of the "nodes method" 
30 according to the invention) . 
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From assumed information about the thickness e of 
the scrutinised layer and about the refraction index n(A) 
and about the extinction coefficient k(A) of the material 
of this layer, we create initial nodes (only a few) and 
5 • an initial thickness e. 

We thus have m nodes and, via interpolation, we can 
find n(A) and k(A) when X is different from the values of 
the wavelengths associated with the nodes. 

From the initial thickness e and these initial 
10 values n(A) and k(A), we determine the theoretical 

spectrum *P by using ellipsometric and ref lectometric 
calculations . 

Furthermore, by means of ellipsometric and 
ref lectometric devices and a spectrometer, we obtain 
is Si, 2 (6,A) and R(A) and we deduce the measurements marked 
W (for measuring conditions 9 and A). 

We then compare W and W by using an error function 
Er and we optimise the refraction index value and the 
extinction coefficient value at the different nodes as 
20 well as the thickness value, by trying to minimise 

Er (Y, W ) . 

Once these values have been optimised and if the 

accuracy on the spectrum n(A), the spectrum k(A) and the 
thickness e is insufficient, we introduce new nodes, we 
25 make the thickness e vary and we recommence the 

determination of V , the comparison of ¥ and W and the 
optimisation that are aforementioned, etc. 

We break the thus defined loop when the accuracy on 
each of the spectrums n(A) and k(A) and on the thickness 
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8 is considered to be sufficient (satisfactory adjustment 
of W and W ) . 

The spectrums n (A) , k(A) and the thickness e are 
thus characterised. 
5 In a non-restrictive illustration we find a layer 

thickness e equal to 212. 3nm. 

Figure 1 diagrammatically shows the scrutinised 
layer 2, created on a substrate 4. We can see the 
ellipsometric device 5,6, the ref lectometric device 8 and 
10 the spectrometer 10. Furthermore, we can see electronic 
processing means 12, comprising a computer and allowing 
to characterise n(A), k(A) and e according to the 
information supplied by the spectrometer 10 and according 
to the method of the invention. 
15 These means 12 are equipped with display means 14 

which allow, in particular, to display the curve of the 
variation of n as a function of A and the curve of the 
variation of k as a function of A. 

We will now consider the illustration in greater 
20 detail. 

Phase 1 

The method of this example firstly comprises an 
initialisation stage . 

The algorithm starts with a limited number of nodes, 
25 more precisely at least one node. We can therefore start 
with a single node, by imposing a refraction index and an 
extinction coefficient which remain constant when the 
wavelength varies. 

We chose node positions so as to, from this node 
30 family, deduce the entire spectrum via interpolation. The 
layer model therefore has 3 parameters or more, as the 
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thickness is also a variable to be determined. The index 
table over the entire spectrum is therefore deduced from 
the nodes via interpolation. 

This is the case when the thickness of the layers is 
5 about the dimension of the wavelength (thin layers) . But 
when the thickness of the layer is greater than the 
coherence length of the light source, the thickness 
barely enters into the calculation of the layer response, 
and is therefore no longer a variable of the problem. 

10 For example, in the case of optical disks (CDROM) , 

on which a very thick layer is deposited (a thickness of 
about 1 millimetre) , the reflection coefficient of such a 
layer is no longer a function of the thickness of the 
layer but purely of the refraction index of the latter, 

15 the coherence length of the incident light beam being 
less than the thickness of the layer. In this precise 
case, the coherence length of the incident beam is 
determined by the roughness of the layers. 

We choose, for example, to place the first two nodes 

20 at the ends X min and X max of the spectrum. The complex 
index values at these ends are chosen according to the 
type of material under scrutiny. For example, on an 
ellipsometric spectrum between 300nm and 800nm of a thin 
layer of photoresist, we take n (300nm) =n (800nm) =1 . 5 and 

25 k (3 00nm) =k (800nm) =0 . 

When the spectrum is only characterised by two 
nodes, the index between the ends is determined via 
linear interpolation. In the case under consideration, we 
therefore have n(X)=1.5 and k(X)=0 for X belonging to 

30 [300nm, 800nm] . 
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From three nodes up, we rather opt for cubic 
interpolation in order to obtain softer forms of index 
law compared to those obtained via linear interpolation. 

The initial thickness is, itself, chosen as close as 
5 possible to the real thickness. 

Phase 2 

We then proceed with an optimal determination of the 
values of the refraction index and of the extinction 
coefficient on the nodes and of the thickness value. 

10 To accomplish this, the spectrums W (A) are 

calculated by means of the used layer model, resulting 
from the choice of nodes, the interpolation law and the 
thickness of the layer. 

The physical model used to calculate ¥ of course 
15 depends on the applied method for measuring, i.e. in 
particular on the incidence angle of the light, on the 
applied spectrum and on the model of thin layers or of 
thick layers if appropriate (see for example the model of 
stacked layers in document 3) ) . 
20 The spectrums W being constituted of a set of 

measurements of diverse nature (for example ellipsometric 
and ref lectometric measurements) , we use a ref lectometric 
(respectively ellipsometric) calculation method for 
ref lectometric (respectively ellipsometric) measurements . 
25 The ref lectometric and ellipsometric measurements 

are combined by means of an error function Er (¥, Y ) 
which is for example of the type defined by the equation 
(1) . 

Thanks to a search function, we minimise the 

30 difference between W (A) and ¥ (A) , by varying the value 
of the refraction index and the value of the extinction 
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coefficient at the position of each of the nodes as well 
as the layer thickness (if this thickness is an important 

factor in the calculation of W ) . When the difference is 

minimal, that being when Er (W , W ) is minimal, this 
5 signifies that the ref lectometric and ellipsometric 
measurements coincide to their best (for a given number 
of nodes) . 

At this stage we obtain, for a known number of nodes 
and a known spectral position for each of these nodes, 
10 the layer model (refraction index, extinction coefficient 
and thickness) which best corresponds to the real layer. 

The validity of the found model is better ensured as 
the number of measurements increases. To obtain a large 
number of measurements we can for example use several 
15 incidence angles of the light 6i, 1 < i < 1, for the 
ellipsometric spectrums, carry out a ref lectometric 
measurement and complementary goniometric measurements. 
Phase 3 

Then, we increase the number of nodes. 
20 We add a finite number of nodes. In a first 

embodiment, the added nodes are positioned so as to best 
represent the spectrums n (A) and k(A). By way of example, 
we place these additional nodes at spots where the 

difference between W and W is maximum or at spots where 
25 the nodes are the furthest apart. We then return to phase 

2 as long as the accuracy on each of the spectrums n(A) 

and k(A) and on the thickness e is insufficient, i.e. not 

equal to a pretedermined accuracy. 

In a second embodiment, we insert new nodes between 
3 0 two nodes of the set of previously chosen nodes, these 

new nodes being evenly distributed over the spectrum. 
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All the same it is appropriate to heed the 
following. 

When we increase the number of nodes, the position 
of the old nodes does not need to be conserved. For 
5 example, assume that we evenly sample a spectrum from 
400nm to 800nm with 3 nodes. These nodes are therefore 
respectively located at 400nm, 600nm and 800nm. When we 
increase to 6 nodes, the three additional nodes cannot be 
placed so that the spectrum is evenly sampled if we 
10 conserve the position of the old nodes. The position of 
the 6 nodes can be defined, is we wish to have even 
sampling, by the values 400, 480, 560, 640, 720 and 
800nm. The old middle node at 600nm therefore disappears. 
To calculate the index value at these positions from the 
15 old nodes we use interpolation. 

In the following we provide two common examples of 
application of the invention. These examples implement 
two types of materials which apply different laws. 

From ellipsometric and ref lectometric measurements 
20 we propose to find the physical laws applied by these 
materials . 

We proceed in the following manner. 

A fictitious material is created, this material 
applying a known theoretical law (a Cauchy law or a 
25 harmonic oscillator law) , with parameters that we 
arbitrarily set. The variations of the complex refraction 
index as a function of the wavelength are thus perfectly 
known. 

Furthermore, we impose a thickness of 200.00nm for 
30 the material, on a silicon substrate, this latter also 
being well known. 
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Fictitious measurements (ellipsometric measurements , 
reflectometric measurements) are calculated, then noise 
induced so as to introduce an apparatus defect. 

It all takes place as if we had used real 
5 measurements taken on the material . But contrary to 
reality we fully know the complex refraction index as we 
set it just as we set the thickness of the layer of 
material . 

Here we '"blindly" test the method, that meaning we 
10 start with an erroneous thickness (22 0nm) and false 
complex refraction indexes, as they are supposed to be 
unknown . 

We apply the method according to the invention then 
we compare the complex refraction index found with the 
15 theoretical complex refraction index. We precisely find 
the same laws, as well as the same layer thickness. 

As a first example, take a material whose complex 
refraction index applies to a Cauchy law such as: 

/ x 300 2 300 4 

n(>l)=1.5 + 0.1 — ^- + 0.1 — -r- 

20 k(A)=0 

This index law is typical of photoresists (in the 
spectral range from 300nm to 800nm) . 

In order to find, using the nodes method (i.e. the 
method according to the invention) , the aforementioned 
25 index law, we carry out two measurements, namely an 
ellipsometric measurement at a 70° angle and a 
reflectometric measurement . 

The processing conditions are' as follows: 
- the processed spectrum lies between 3 00nm and 
3 0 8 0 0nm; 
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- at the start, the nodes are in the positions 
(300nm, 1.6) and (800nm, 1.6), this meaning that the 
index is considered as linearly varying between 3 00nm and 
800nm # and that its value is constant (equal to 1.6); 

5 - the number of nodes is repeatedly increased 

according to the sequence 2-»4— »6; 

- during the node increasing procedure, the 
wavelength position of each node is calculated so that 
the sampling by 1/A is regular (X: wavelength) , the 

10 density of points thus being increased around the low 
wavelengths ; 

- the interpolation law is a cubic law, when the 
number of nodes is greater than 2, otherwise it is 
linear; and 

15 - the applied minimisation algorithm is a Simplex 

type algorithm. 

Figure 2 allows to compare the refraction index 
corresponding to the fictitious material which perfectly 
applies the Cauchy law (curve I) with the refraction 

20 index which we find via the nodes method (curve II) , 
using 6 nodes (represented by circles in figure 2) . We 
have used a W composed of an ellipsometric measurement 
{Si (X) , S 2 (X) } at 70° and a ref lectometric measurement 
R(X). 

25 We finally find a thickness of 199.8384nm. 

As for a second example, take a material whose 
complex refraction index applies a law with two harmonic 
oscillators, such as: 



30 



[n(E)+jk(E)] 2 =l+ £ 



i^E + Ej + jGj E-Ej + jGj 
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with 



j 2 



= -1 and E 



1240/A (A in nm) 



A 1= 0 .25x1 . 5 2 



A 2 =0.25xl.5 2 



E!=1240/400 



E 2 =1240/300 



Gi=0.3 



G 2 =0.3 



10 



15 



20 



25 



In this second example, the nodes method is applied 
to a set of ellipsometric measurements carried out 
between 250nm and 800nm, at 75°, 70°, 60° and 45° , along 
with a ref lectometric measurement. The real thickness of 
the material being 200nm / we find a thickness of 200.25nm 
with the nodes method. The adjustment to the index law 
considered in this second example is very good, as 
figures 3A and 3B demonstrate. 

These figures 3A and 3B respectively illustrate the 
reconstructions of the curves n(A) and k(A) of the 
material by means of the nodes method. The reconstruction 
is carried out using four ellipsometric spectrums and a 
ref lectometric spectrum. The real absorption peaks are 
particularly well represented by the curve obtained via 
cubic interpolation between the nodes (represented by the 
circles in figures 3A and 3B) . 

In figure 3A, the curve I (respectively II) 
corresponds to a refraction index n which perfectly 
applies the chosen law (respectively to a refraction 
index n found via the nodes method) . 

In figure 3B, the curve I (respectively II) 
corresponds to an extinction coefficient k which 
perfectly applies the chosen law (respectively to an 
extinction coefficient k found via the nodes method) . 

Examples of the invention have just been described 
above. More generally speaking, we note that in the 



B 14206.3 PV 



26 



invention we consider a set X of values, with X= 
{n 1/ n 2# .« / n i/ ..Ji m ,k 1/ k2,...ki / ...k m/ e} , where 

ni is the value of the refraction index (real) 
corresponding to Ai, ie{l...m}, m being the number of 
5 nodes 

ki is the value of the absorption coefficient at the 
node corresponding to Ai, ie{l...m}, 

e is the thickness of the scrutinised layer. 

In this case, the operation for minimizing the error 

10 Er (W, W ) means finding the set or "vector" X such that 
Er is minimal. 

When we do not assert any specific constraints, the 
minimising is a minimising with 2xm+l parameters. We can 
of course introduce constraints in order to reduce the 
15 number of variables. 

In particular, if we know that the material is non- 
absorbent, we impose ki=0 for all i of {l...m} and X 
becomes : X= { n x , n 2 — , n± , ...n m/ e } . 

If through a complementary measurement (for example 
20 a goniometric measurement or a non-optical direct 
measurement) we know with sufficient accuracy the 
thickness of the scrutinised layer, the thickness e is no 
longer a variable and we obtain: 

X= { rii , n 2 , ... , ni , ...n m , ki , k 2 , ...ki , ... , k m } . 
25 Of course, the two previous options can be combined. 

In the following we explain an embodiment of the 
invention in a more general form than the previous 
examples . 

Given that a (A) is a function of the wavelength A 
30 of the light used in the measurements, we can for example 
choose : 



B 14206.3 PV 



27 



a(A)=A (see figures 3A and 3B where the spectrum is 
evenly sampled in X) 

a(A)=l/A (see figure 2 where the spectrum is evenly 
samples in 1/A) 

5 a(A)=hc/A where h is the Planck's constant and c the 

speed of light in the vacuum, oc{\) then being homogenous 
with an energy. 

Given that A is the interval of the spectrum of 
measurement, B is the spectral interval described by the 
10 "nodes" and C is the interval of interest. 

The interval C is included in the interval B or 
equal to this interval B. Likewise, the interval B is 
included in the interval A or equal to this interval A. 

We specify that each of the intervals A, B and C is 
15 of [a m , a M ] type where a m is less than a M and there are 
two wavelengths A k and A± such that a m =a (A k ) and a M =a(A!). 
By way of non-restrictive illustration, figure 4 shows an 
example of the intervals A, B and C and of the variation 

curves of ? and W as a function of a (A), W in fact being 
20 a function of n*(a(A) ) - 

The circles N represent the nodes. 

n* is the complex index which is expressed here as a 
function of a (A) and whose real and imaginary parts are 
respectively marked n(a(A)) and k(a(A)). 

25 We also see, by way of non-restrictive example, an 

example of the variations curve of n (respectively k) as 
a function of a (A), passing through the points having co- 
ordinates (ni, aj.) (respectively (k if oc±) , where ai=a(Ai), 
ni=n(ai), ki=k(ai), l<i<m (m being a natural number other 

30 than zero) . 
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In figure 4 we note the relation between the nodes 
N, the points (n±, oc±) and the points (k if oci) . 

In the embodiment under consideration, we use an 
algorithm comprising the following stages: 
5 1 . we carry out the measurements W over the 

interval A and we choose the calculation methods 
associated with the measurements (ellipsometric or 
ref lectometric calculations) ; 

2. we choose m numbers ot± (constituting m initial 
10 values of the function a), i belonging to {l,..jn}, with 

m>l, and {cx±} <= A (the oc± corresponding to the "nodes"); 
when m>l, B is defined as the set of points a such that 
min(ai) <a< max(ai); when m=l, we have B=A; 

3. we choose m initial values of complex index n* 
15 at the m points ai, i belonging to {1,..jti}; 

4. if m^l we choose an interpolation law which 
allows to calculate the refraction index n(a) over the 
interval B from the points (a i# n ± ) , i belonging to 
{l,...,m}; if m=l, then n(a)=ni(ai) over the entire interval 

20 B; 

5. we choose M variable parameters with M<2m+1; 
these parameters can be for example : 

— the real parts of the refraction indexes at 
points ct±, i belonging to {l,...,m}, or 

25 - the imaginary parts of the refraction indexes 

at these same points, or 

- the thickness of the material for which we are 
searching the refraction index; 

6. we choose an error function Er , W) which 
30 characterises the difference between a measured spectrum 
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and a theoretical spectrum; generally, the error is 
measured over the interval C; 

7. using a minimising function of Er (W , W) with M 
parameters, we perform the following series of stages: 

5 a) by applying the interpolation function of 

(aj^ni) over B, we deduce n(a) with a belonging to B; 

b) by using n(a) and the thickness e and methods 
for calculating spectrums, we calculate the theoretical 

spectrum ¥ (n(a),e); 

10 c) we compare W and W by using Er=Er (¥, W ) ; if Er 

is sufficiently small (i.e. if Er is less than a 
predetermined value e) , or if Er is minimal, we go to 
stage e) , otherwise we go to stage d) ; 

d) we make the M variable parameters vary so as to 
15 tend to the minimum of Er and we return to stage a) ; 

e) if Er is less than e we then obtain a set of M 

parameters such that Er {W, W (n (a, M) , e ) ) is minimal and 
the refraction index calculation is completed: this index 
is taken equal to the last one obtained; and if Er is 
20 greater than or equal to e we go to stage 8) . 

8. We increase m and we return to stage 2) . 

It is appropriate to note that the invention can not 
only use wavelength (A) samplings but also frequency (c/A 
or simply 1/A) samplings, energy (hc/A) samplings and, 

25 more generally, samplings with parameters which are 
functions of the wavelength. 

Furthermore, it is appropriate to note that an 
aforementioned essential stage of the algorithm (stage 8) 
is not restricted to adding a set of nodes to the already 

30 existing nodes: it encompasses the more general case 
where the number of nodes increases. This means, in a 
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specific embodiment of the invention, that after 
minimising with 3 nodes, if we want to increase to 6 
nodes in all, the position of the 3 old nodes is "erased" 
so as,' for example, to have a constant density of nodes 
5 over the spectrum. In practice, this is the best option. 
The information about the position of the old nodes is 
not lost as the index value at the old nodes is used to 
calculate the values of the 6 new nodes (in fact 3 new 
nodes plus 3 old nodes) . 

io Thus, according to a specific embodiment of the 

method according to the invention, we can increase the 
number of initial values of the function a by adding one 
or several values to the extant initial values; however, 
according to a preferred embodiment, we can increase the 

15 number of initial values of the function a by replacing 
the extant initial values with new initial values whose 
number is greater than the number of extant initial 
values . 

The invention in not restricted to the 
20 characterisation of thin layers. It also applies to the 
characterisation of thick layers. 

Furthermore, the invention is not restricted to the 
characterisation of a single layer, created on a 
substrate. It also applies to the characterisation of 
25 two, or more than two, layers created on a substrate. 
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